Spin-dependent diffraction at ferromagnetic/spin spiral interface 

A. Manchon 
SPIN TEC, URA 2512 CEA/CNRS, 
CEA/ 'Grenoble, 38054 Grenoble Cedex 9, France 

N. Ryzhanova 
SPINTEC, URA 2512 CEA/CNRS, CEA/ Grenoble, 
38054 Grenoble Cedex 9, France and 
Department of Physics, M. V. Lomonosov Moscow State University, 119899 Moscow, Russia 

A. Vedyayev 

SPINTEC, URA 2512 CEA/CNRS, CEA/ Grenoble, 
38054 Grenoble Cedex 9, France and 
Department of Physics, M. V. Lomonosov Moscow State University, 119899 Moscow, Russia 

B. Dieny 
SPINTEC, URA 2512 CEA/CNRS, 
CEA /Grenoble, 38054 Grenoble Cedex 9, France 
(Dated: February 2, 2008) 

Abstract 

Spin-dependent transport is investigated in ballistic regime through the interface between a fer- 
romagnet and a spin spiral. We show that spin-dependent interferences lead to a new type of 
diffraction called "spin-diffraction". It is shown that this spin-diffraction leads to local spin and 
electrical currents along the interface. This study also shows that in highly non homogeneous mag- 
netic configuration (non adiabatic limit), the contribution of the diffracted electrons is crucial to 
describe spin transport in such structures. 

PACS numbers: 72.25.-b,73.43.Qt,75.60.Ch,75.50.Ee 
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The recent observations of current-induced domain wall motion [1] and the investigation 
of this phenomenon by micromagnetic simulations [2] have underlined the question of spin 
transport in non homogeneous systems. Although proposed very early [3], current-induced 
domain wall motion (DWM) has attracted much attention because of its great application 
potential but also because of the deep fundamental question of the role of electron spin 
motion in temporally and spacially varying magnetic structures. Zhang et al. [4] proposed 
an analytical formula for spin torque in slowly varying structures (see eq. 9 in Ref. 4) : 
dM dM 

T = a— + bM x — - Cl M x[M x (j e .V)M] - c 2 M x (j e .V)M (1) 

where M is the magnetization unit vector, a is the renormalization factor of the gyromagnetic 
ratio and b is the renormalization factor of the damping parameter in the Landau-Lifshitz- 
Gilbert equation. The last two terms represent the spin transfer torque in spacially varying 
magnetic structure. The prefactors c\ and c 2 are respectively proportionnal to the adiabatic 
and non adiabatic contribution. For smooth enough non-homogeneities of the magnetic 
structure, the adiabatic approximation is usually assumed : the electron spin follows the 
local magnetization producing a small torque proportional to the spatial derivative of the 
magnetization [5, 6]. In this case, it is usually accepted that the non adiabatic term c 2 is 
small but cannot be neglected in domain wall experiments [7, 8]. Furthermore, in spacially 
non homogeneous systems, C\ and c 2 are non local coefficients [5, 6, 9]: at each point of 
the structure, one has to consider the contribution of all the electrons flowing through the 
structure. 

In this article, we propose to study the non adiabatic regime in an "academic" system. 
We consider two adjacent layers: the left one (z<0, F) with a homogenous magnetization 
P = Pz and the right one (z>0, SS - for Spin Spiral) with a 2D helical magnetization 
contained in the (x,z) plane M = M(sin9(x)x + cos 9(x)z), where 9(x) = 9 + Qx. The 
interface lies in the (x,y) plane and z is perpendicular to the interface. We consider that 
the regions are semi-infinite and respectively connected to a ferromagnetic and a spin spiral 
reservoir. The bias voltage V is applied across the interface. 

Such helical spin structures exist in some compounds such as MnSi [10], oxide materials 
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such as SrFeOs, NaCuO [11], rare-earth based compounds [12] or 7-iron [13]. This helical 
structure can be also a simplified picture of narrow stripe domains with domain wall width 
comparable to domain width. 

In this case, the adiabatic approximation is no more valid because an electron moving 
through the interface keeps the memory of its spin state for some distance. Spin-polarized 
electrons moving from F into SS undergo "spin diffraction" as represented on Fig. 1: an im- 
pinging electron with an in-plane incident wavevector k gives rise to transmitted (reflected) 
waves with in-plane wavevectors n + nQ/2 (n+nQ), 2ttQ~ 1 being the wavelength of the spin 
spiral. From the interference of all these waves, one can expect original electrical properties 
like non-zero torque acting on the magnetization of SS, local longitudinal spin current and 
even charge current (Hall effect and spin Hall effect [14]) along the F/SS interface. 

To model the spin-dependent transport, we use Keldysh out-of-equilibrium technique[15, 
16] which expresses the lesser Keldysh Green functions G~^(rr') as a function of the basis 
of wavefunctions ^f^ rr7 (r) for an electron moving from the left (right) to the right (left) 
reservoir: 

G;+(rr') = f^Ur^Ur) + fMKA^Uv) (2) 

where fi( r )(^i(r)) are the Fermi distribution functions in the left and right electrodes and 
fj,i( r ) are the chemical potentials in these electrodes so that V = (fj,i — fi r )/e. The electrical 
current density J e , spin density m and torque T (exerted on SS) are given by the usual local 
definitions: 

je = 3 [(Vr , _ V r )(Gf+(rr') + <V(rr'))] (3) 

m x + im y =< o + >= 2G^ + (rr') (4) 
m z =< a z >= Gf T + (rr') - G u + (rr') (5) 
T = x m (6) 

hfjL B 

where 7 is the gyromagnetic ratio, \ib is the Bohr magneton and J s d is the s — d exchange 
coupling. The torque T possesses two components: the usual spin transfer torque term 
(STT, adiabatic torque propotionnal to m y and lying in the (x,z) plane), and the field-like 
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term (IEC, non adiabatic torque perpendicular to the (x, z) plane). The Hamiltonian of the 
system is then: 

2 

H=^- + U + J sd {a.S) (7) 

where U is the potential profile, cr is the vector of Pauli matrices and S is the magnetization 
of the layer (S = M in SS and S = P in F). The wavefunctions are the solutions of 
the Schrodinger H^/ = E^/. To solve these equation, we used the procedure developped 
by Calvo[17] in a spin spiral. The boundary conditions of these wavefunctions and their 
derivatives lead to recurrent relations between the coefficients of these wavefunctions [18]. 
For example, the wavefunctions of an initially majority electron originating from the left 
reservoir with in-plane wavevector k are: 



*S8 = E[C3,. +0 ,„- 4 »(COB - i0„ +e( „_., sin ?|V** + «"-*>>* 

n 

+ ^(n-fcjC W"*> COS ^ + i sin ^) e «^(-l))-] e «(«H-(«+«(«-i»*) ( 8 ) 

*S8 = E[^(n- § ))(sin °M + cos e M )e ik ^-»* 

+ ^(n-f ))H cos ^ + K+Q(n _ f) sin ?M )e *^(-l))* ]e ^(^(»-i)W (9) 

where 03(4) are complex diffraction coefficients and ^3(4) are wavevectors for majority (mi- 
nority) spin projection in SS. This procedure will be developped in a forthcoming article 
[18]. 

For the numerical simulations, we took parameters corresponding to spin transport in 
Co: the Fermi wavevectors for majority and minority spins are respectively k' F = 1.1 A -1 , 
4 = 0.6 A -1 ; the inverse wavelength of the spin spiral is Q 1 = (2tt) 1 A 1 (highly non 
homogeneous magnetic system). The Fermi wavevectors of SS and F are the same and 
we consider the linear approximation : for a small enough bias voltage, only the electrons 
originating from the left reservoir with an energy located between n\ and fi r significantly 
contribute to the charge and spin transport. We set fii — fi r — 38 meV. 

Figure 2 displays the longitudinal charge current J% in the (x, z) plane. Close to the 
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interface, oscillates with SS magnetization, but rapidly decreases and vanishes to zero 
within 5 A. We observe the same behaviour for the perpendicular current J e z except it reaches 
an averaged value within 5 A. We observe the same features for STT and IEC (not shown 
here). 

This rapid decay is due to interferences between diffracted electron waves. Close to the 
interface, the averaging effect is smaller than in the bulk SS, so the oscillation amplitude of 
J%{z) and STT, IEC is always higher near the interface than in the bulk. 

These damped oscillations strongly depend on the SS wavelength 27r<5 _1 [18]. Fig. 3 
displays the ^-dependence of the longitudinal current J% [Fig. 3(a)], STT [Fig. 3(b)] and 
IEC [Fig. 3(c)] for different values of SS wavelenght Q _1 . 

When Q increases, the oscillation amplitude of J| is reduced so that vanishes more 
rapidly to zero: the increase of the non adiabaticity induces a more important averaging 
effect due to interference between multiple diffracted waves. STT and IEC have an opposite 
behaviour when varying Q. The oscillation amplitude of STT decreases when Q increases 
(similarly to JJ), whereas the amplitude of IEC increases. This illustrates the different 
nature of STT and IEC: STT is the adiabatic torque (vanishes in highly non homogeneous 
structure) and IEC is the non adiabatic torque (vanishes in adiabatic systems). 

This study demonstrates that spin diffraction gives rise to complex characteristic in spin 
torque and electrical currents and is of seminal importance in non adiabatic magnetic sys- 
tems. 

This work was partially supported within the European MRTN SPINSWITCH CT-2006- 
035327 and the Russian Fundings for Basic Research 07-02-00918-a. 
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FIGURE CAPTIONS 




Figure 1: Cartoon of the bilayered structure. The left semi-infinite layer is a ferromagnet with 
a homogeneous magnetization and the right semi-infinite layer is a spin spiral with wavelength 
2irQ~ 1 . The spin-polarized electrons undergo spin-diffraction at the interface. 




Figure 2: Map of the longitudinal electrical current density in the (x, z) plane. 
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Figure 3: Longitudinal current (a), usual spin transfer torque (b) and current-induced interlayer 
exchange coupling (c) as a function of z for different Q (see inset) and calculated at x = tt/2Q. 
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